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I , Abstract 

Ph. 
(D . 
^ • We present a renormalized computational framework for the evolution of a 

self-interacting scalar field (infiaton) and its quantum fluctuations in an FRW 
K> ! background geometry. We include a coupling of the field to the Ricci scalar 

5j I with a general coupling parameter ^. We take into account the classical and 

quantum back reactions, i.e., we consider the the dynamical evolution of the 
cosmic scale factor. We perform, in the one-loop and in the large- A^ approx- 
imation, the renormalization of the equation of motion for the infiaton field, 
and of its energy momentum tensor. Our formalism is based on a perturba- 
tive expansion for the mode functions, and uses dimensional regularization. 
The renormalization procedure is manifestly covariant and the counter terms 
are independent of the initial state. Some shortcomings in the renormaliza- 
tion of the energy-momentum tensor in an earlier publication are corrected. 
We avoid the occurence of initial singularities by constructing a suitable class 
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of initial states. The formalism is implemented numerically and we present 
some results for the evolution in the post-inflationary preheating era. 



1 Introduction 

Nonequilibriuin processes in cosmology have recently been considered by var- 
ious authors. The main interest has been centered around the possible infla- 
tionary period of the universe (see e.g. [|l], ^, ^) and the subsequent reheating 
0, H . It has already been found, by considering parametric resonance ^ as- 
sociated with oscillations of the inflaton field and by exact computations 
including back reaction both in Minkowski space |0, H and in an expand- 
ing universe |^, 0, |1T], |12[, that the time-dependent inflaton field produces 



particles or classical fluctuations preferentially at low momenta and not in a 
distribution corresponding to thermal equilibrium. This process of particle 
production has therefore been termed preheating |^ |13]. 

The equations of motion for nonequilibrium systems have been presented 



by various authors [|I4[ ^ using the CTP formalism introduced by Schwinger 



jT6| and Keldysh ||T^ . Their application to inflation within a conformally flat 



Friedmann-Robertson- Walker (FRW) universe has been initiated by Ring- 



wald [|T8|; they have recently been implemented numerically in ^, |T^, |T^ 
and |12[. Similar computations have been performed in configuration space 
for the case that the fluctuations are treated as classical ones [^. Again 
fluctuations with rather low momenta are strongly excited, thus justifying 
the classical approximation. Apart form such exact numerical computations 
there exist also various analyses based on analytical approximations to the 
solution of the mode equations [|lT], ^, 0, |2^, ^ . 

We have recently considered [^ nonequilibrium dynamics of a scalar 



inflaton field and its quantum fluctuations in FRW cosmology. This work 
was based on a method for renormalized numerical computations in quantum 
field theory, introduced for nonequilibrium dynamics in Minkowski space in 
||26|| . In this method the leading, divergent parts of the fluctuation integrals 
appearing in the inflaton equation of motion and in the energy-momentum 
tensor are separated from the numerical computation of finite parts. This 
allows for a free choice of regularization, which is performed analytically. In 
p5| we chose dimensional regularization. The renormalization then can be 
performed in the usual way, with the standard counter terms of equilibrium 
quantum field theory, manifestly covariant f\ and independent of the initial 
conditions. In applying this method to a scalar field in FRW space-time we 
found that renormalization introduces singularities in the time variable at 

''This is to be understood in the restricted sense of special relativistic covariance. 



the initial time to, taken to be in the following. While the equation of 
motion for the inflaton field is finite as t ^ 0, the energy- momentum tensor 
is not. So the Friedmann equations have an initial singularity. We, therefore, 
were not able to numerically implement our formalism. 

Depending on the parameters and the model under consideration, these 
singularities can be more or less pronounced. It is well known that the 
main excitation of fluctuations appears in a low momentum resonance band, 
so the ultraviolet divergences, the precise handling of renormalization and 
thereby the initial singularities may be considered relatively unimportant. 
Certainly the main features of inflaton dynamics, as presented in |]I0[ will not 



be changed if renormalization is handled in a more meticulous way. However, 
in order to be on safe grounds, one should be sure that the renormalization 
aspects can be handled properly and consistently. 

In |T^ renormalization was performed by adiabatic regularization, i.e.. 



by subtracting the leading adiabatic orders of the fluctuation terms. The 
problem of initial singularities there was avoided by starting the system with 
a smooth transition, for the inflaton mean field and its fluctuations, from 
an initial phase in which the back reaction of the fluctuations is switched 
off. In this initial phase the time dependence of cosmic scale parameter is 
fixed, and the evolution of the fluctuations is handled analytically. The initial 
singularities are indeed related to a noncontinuos behaviour of the effective 
mass of the scalar field fluctuations, a mass whose variation is determined 
by the background metric and the inflaton mean field. In |^ (see also 
||28|| ) we found a different solution of the problem of initial singularities, a 
Bogoliubov transformation of the initial state. Using our mode expansion 



and our construction of the initial state, it was shown in |2^ that this initial 
state corresponds, in the terminology of the adiabatic expansion, to a vacuum 
state of adiabatic order 4. 

Having solved the problem of initial singularities we are now able to im- 
plement our formalism numerically. Besides these numerical computations 
with nonsingular initial conditions we will present here some further formal 
developments. 

Part of these is made necessary by fact that in our previous work we 
omitted several finite terms in the renormalized energy-momentum tensor, as 



criticized in [^ . In applying dimensional regularization we did not take into 
account the dimensional continuation of the conformal rescaling of the fields, 
and of the dimensional continuation of the various basic tensors. Further- 
more, the determination of the counter terms was based on the consideration 



of the conformally flat FRW metric. It is not possible, then, to fix these 
counter terms in an unique way. As a consequence the anomaly of the stress 
tensor did not appear at all. This is, therefore, not due to a shortcoming of 
our formalism, but to the very special nature of the problem. 

A further extension of our previous work is the consideration of the large- 
A^ limit in the 0{N) a model. This requires a reconsideration of the renor- 
malization procedure, leading to modified counter terms. 

The analysis presented here is performed for arbitrary values of the con- 
formal coupling ^, which appears in various renormalization constants and is 
itself renormalized. 

We will present the basic relations of FRW cosmology in section 2. The 
nonequilibrium dynamics of a scalar mean field and its fluctuations in this 
geometry is introduced in section 3, the associated energy-momentum tensor 
in section 4. In section 5 we describe our perturbative expansion of the mode 
functions and derive some expressions occuring in the various fluctuation 
integrals. Renormalization of the one-loop equations of motion and of the 
energy-momentum tensor is considered in sections 6 and 7, respectively. In 
section 8 we construct the Bogoliubov transformation by which the initial 
singularities are removed. In section 9 we extend our formalism to the 0{N) 
a model. We present and discuss some results of our numerical computations 
in section 10. 

2 FRW cosmology 

We consider the Friedmann-Robertson- Walker metric with curvature param- 
eter /c = 0, i. e. a spatially isotropic and fiat space-time. The line-element 
is given in this case by 

ds^ = dt^ - a^{t)dx^ . (2.1) 

The time evolution of the cosmic scale factor a{t) is governed by Einstein's 
field equation 

(1 + 6Z)G^, + <5a«if ^, + 5(3^^^ H^, + s^ h ^, + 5Ag ^, = ~k{T^,) (2.2) 

with K = SttG. 

The Einstein curvature tensor G^^ is given by 
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(2.4) 
(2.5) 



The Ricci tensor and the Ricci scalar are defined as 

R = g^'^R^u , 
where 

^a/S-r = '^■y^al3 ~ '^a^-y/3^ -jamais ~ ^aa^-r/S ■ (2-6) 

The tensors ^^^H^^, ^"^^H^^^ and H^i, arise from the variation of terms pro- 
portional to R^.R'^^Rap, and R°'^'^^Ra/3-ys in the Hilbert-Einstein action. 
Their precise definitions, valid for general dimension n, and the subsequent 
identities for n = 4 are given in [RO]: 
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- --g^uR'^'^^'RaMS + 2R^^^^-Rf^ - AUR^, + 2R,^, 
-m^^RZ + m'^^Ra^pu. (2.9) 

In the case n = 4 the generalized Gauss-Bonnet theorem states that 



Jd^x^/^ (r^p^sR'"^^^ + R^- ^RapR'"^) 
is a topogical invariant. It then follows that 
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flU 



«i7„. + 4(2)if 



flU 



fif ; 



n 



Furthermore, in conformally flat space-time as considered here. 
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(2.10) 



(2.11) 



(2.12) 



These terms, as well as the terms SZG^i, and SAg^^i^ are introduced for the 
purpose of renormalization. They are going to absorb the divergences arising 
in the energy-momentum tensor. It will be more convenient, later on, to add 
analogous terms to the energy-momentum tensor, with coefficients 5A = 
SA/k, 6a = 6a/ n etc.. Of course the cosmological constant and higher 
curvature terms may be present in the bare action already, if required by 
observation. Dividing (|2.2| ) by 1 + 6Z one sees that 6Z may also be considered 
as renormalizing Newton's constant. 

With the FRW metric ( p.l|) the Einstein field equations reduce to equa- 
tions for the time-time component and for the trace of G^^^, the Friedmann 
equations 



(1 + 6Z)Gu + a ^'^Hu + 6P ^'^Hu + 6^Hu + 6K 
;i + 6Z)G^^ + a ^^^Hji + 6(3 (^^i//^ + 6-fHll + n6K 






(2.13) 
(2.14) 



We now compute the various terms for the n dimensional line element 
( |2.1| ). For the Christoffel symbols for an n dimensional fiat FRW universe 
one finds 



J- j^j CiCi O-ij J 
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(2.15) 



The nonvanishing components of the Riemann tensor are 



-^ itk ~ ^'^ ^ik 
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(2.16) 
(2.17) 



for those of the Ricci tensor one finds 

a 



Ru = {n-1)- 



R 



a 



*j 



—da — (n — 2)d^ 



JlJ ; 



(2.18) 



this leads to the Ricci curvature scalar 



R = 2(n- 1)- + (n- l)(n - 2) (- 
a \a 



(2.19) 



Expressed in terms of Bubble's constant 



Hit) 



d(t) 
a(ty 



(2.20) 



it takes the form 

R={n-l)(2H + nH^) . (2.21) 

The time-time components and the trace of the tensors ^'^''^H^^ are given by 

^^'^Htt = -6HR + -R'^ -6H^R (2.22) 

+ (n - 4) [-2HR - (n + 1)RH^) , 

(2)iJ tt = -2HR + -R^ - 2H^R + (n - 4) 1--HR (2.23) 

^ ^{n + 2)H^R + -{n-l){n-2fH^\ , 
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//« = -2HR + -R'^ -2H^R (2.24) 



WiJj; = -6/2-18iJi?+(n-4) ('-2/2-2(n + 2)iJ/?-ii?2') (2.25) 
(2)iJM = _2R-6HR+{n-4)(--R--{n + 3)HR (2.26) 

""^^ + ]{n - 2fH^R - -n{n - l){n - 2fH^\ , 



(n-1) 4 



ff^ = -2i? - 6/7i? + (n - 4) -2iJi? 



2(n 



+ {n-2)H^R--n{n-l){n-2)HA . (2.27) 

3 Nonequilibrium equations for the scalar field 

The Lagrangian density of a 0^-theory in curved space-time is given by 

C = v/^{|9^$9^$ - |m2$2 - |i?$2 _ >L$4| ^ (3,1) 

where R{x) is the curvature scalar and ^ the bare dimensionless parameter 
describing the couphng of the bare scalar field to the gravitational back- 
ground. We split the field $ into its expectation value and the quantum 
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fluctuations ip: 
with 



$(x,t) = 0(t)+7/>(x,t), 



(3.2) 



(3.3) 



Trp(t) ' 

where p{t) is the density matrix of the system which satisfies the Liouville 
equation 

.dp{t) 



dt 



[Hit), pit)] 



(3.4) 



The one-loop equation of motion of a scalar field with A0^ interaction has 



been obtained in the FRW universe by Ringwald |18[; we follow closely his 
formulation. The equation of motion for the classical field is 



+ (n - l)H^ + (m^ + iR)(P + \(P^ + ^(t^^) 



-- . (3.5) 

The expectation value of the quantum fluctuations {ip'^) can be expressed as 



(^2^ = -iG(t,x;t,x) (3.6) 

in terms of the non-equilibrium Green function G{t^ x; t', x') which satisfies 



^ + (n - l)iJ^ + a-'{t)V' + m^ + ^R{t) 



+\^\t) 



G{t,x;t',x.') 



a3(t) 



5(t,x;t',x') . 



(3.7) 



The boundary conditions for this Green functions will be given below. Due to 
the presence of the term H{t)d/dt the differential operator on the left hand 
side of this equation is non-hermit ian. It is made hermitian by introducing 
conformal time and appropriate scale factors. Conformal time is defined as 



dt' 



1 







a{f) 



(3.8) 



In conformal time the line-element ( p.l|) reads 



ds' = C{T){dT' - dx' 



(3.9) 



where the conformal factor C{t) is given by 

C{r) = a\r) . 



(3.10) 



We further rescale the scalar field and its quantum fluctuations by introduc- 
ing the dimensionless 'conformal' fields 



V^(r,x) = at-i(t)^(t,x). 
The Green function is rescaled accordingly via 

^(x, r; x', r) = a5-i(t)at-i(t')G'(x, t; x, t') . 
The equation of motion of the classical field </^(t) now becomes 



(3.11) 
(3.12) 



(3.13) 



^"ir) + a'ir) m^ + (^ - ^ i?(r) v^(r) + 



A(a(r)/i)' 
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with 



e. 



n-2 



ip^r) = (3.14) 



(3.15) 



4(n-l) ' 

and where the primes denote derivatives with respect to conformal time. The 
two-point-function G now satisfies 






\/^ + M^(t) 



G{x, t; X, r) = -(5(x, r; x, r) . 



(3.16) 



Here M'^{t) denotes the square of the effective mass term the fluctuation 

field 'ijj{t, x) 



M^(t) = a^ir) 



m2 + (e-a)i?(r) 



+ 



A(a(r)/i)' 



'V\r) . 



(3.17) 



For later discussion it is useful to divide M'^{r) into the usual 4 dimensional 
part and into a part proportional to (n — 4): 



M'(T)=a' 



m'+[^-^]Rir 



\{afiy 



n- A 
12 in -I 



-a^R. (3.18) 



When multiplied by a divergent factor l/(n — 4) the last term yields a finite 
contribution. 



The problem of determining the Green function is now essentially reduced 
to the equivalent problem in Minkowski space. We expand the fluctuation 
field in terms of the mode functions Uk{T) exp(ikx) via 



V^(r,x) 



d^'-^k 



c{k)Uk{r)e'''^ + c^kM{r)e 



-ikx 



(27r)"-i 
The functions Ukir) satisfy the mode equation 



U'^r) + nliT)Uk{ 



with 



Ql{T) = e + M\r) 



We further impose the initial conditions 

with 



nk{o) = Jk^ + Mi. 



(3.19) 

(3.20) 
(3.21) 

(3.22) 
(3.23) 



In the following we will use the short notation ^2^0 = fifc(O). The nonequilib- 
rium Green function Gfc(r, x;t',x') can by expressed in terms of the mode 
functions via 



G'fc(r,x;r',x' 



Jn— 1 



k 



[r - r')U,{r)Ul{r')e 



V-r'^I^^^Cx-x') 



(27r)-i 21]feo L 

+e{r' - r)Ul{r')Uk{r)e-'^^^-^'A .(3.24) 



The expectation value of the fluctuation fields is given, therefore, by the 
fluctuation integral 



T{r) = {i,\T)) = -iG{t,^-t,^) 



d^-'k |f/,(r)|2 



(27r)"-i 2fifco 
The unrenormalized equation of motion of the inflaton fleld reads 



if" + a^ m'^ + {^-^n)R ^ + 



A(a/i)' 
6~~ 



V' + 



\{ajiY 



^J^=0 



(3.25) 



(3.26) 



The regularization of the fluctuation integral and the renormalized form of 
this equation will be discussed below. 
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The field expansion ( 3.19| ), together with the equation of motion for the 
mode functions and the initial conditions, defines a Fock space in which the 
initial quantum state or density matrix can be represented. The ground state 
is the conformal vacuum state ||3^, it has been chosen previously [|18|, |^ 



as the initial state for the nonequilibrium evolution. As mentioned in the 
introduction we have found [^ that for such an initial state the energy- 
momentum tensor becomes singular at the initial time. The construction of 
a suitable initial state, avoiding such singularities, will be given below. 



4 The energy-momenturQ tensor 



In order to formulate Einstein's field equation we have to discuss the energy- 
momentum tensor of the scalar field in curved space time. For a classical 



field it reads |^ 



T, 



flU 






(4.1) 



In the conformally flat FRW metric the energy-momentum tensor is diagonal. 
One obtains for its time-time component and its trace 



rpcl 
''-tt 



rp^ cl 



1^2 ^ 1^2^2 ^ A ^4 _ ^^^^^2 ^ 2(n - l)^H<P4> 



n 



l-- + 2(n-l)e 



^2 + ^ ^^2^2 ^ ^^4 j _ ^^M02 



24 



+2{n - 1)^ '(j)+{n- l)H(p^ (p . 



(4.2) 



We again introduce conformal time and the conformal rescaling of the fields. 
Furthermore, we include the quantum fluctuations of the field y?. The classi- 
cal energy density then takes the form 



t; 



cl 



a 



^|2^-"4'"V.^ 



<^ 



+2{n - 1) (e - in) ( -W' - \{n - 2)i7V' 
a 4 



(4.3) 



*We continue to consider Tu instead of Trr for convenience. 
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The fluctuation energy density is given by 



n 



d'^-^k \U, 



/|2 



1 



" a2-^ J (27r)"-i [ 2a2 ' 2a^ 

+ (n-l)(e-e 



Tr^nJrYlUk 



(4.4) 



a dr^ ^ 1 2 ' 2 n - 1] 



\U, 



We obtain for the classical and fluctuation parts of the trace 

l2 



rpclfj, 



2-e 
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2(n-l)(^-en) 



^-^(n-2)Hcp 



+2{n-l)^^- + n 



1 2 2 , KafJ-Y 4 

-m 09 H --09 

2 ^ 24a2 ^ 
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2(n-l)e 
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(4.6) 



Energy density and pressure are related to the energy-momentum tensor via 



Ttt 



£ 



Tj: = S-{n-l)p. 



(4.7) 



It is straightforward to show, using the equations of motion for the clas- 
sical fleld and for the mode functions ( p. 7] ), that the energy is covariantly 
conserved: 

£'{T)/a{r) + {n- 1)H{t){p{t) + £{t)) = . (4.8) 

5 Perturbative expansion 

In order to prepare the renormalized version of the equations given in the 
previous section we introduce a suitable expansion of the mode functions, 
which was used in PSI, 0, 15^ for the inflaton fleld coupled to itself, to 



gauge bosons, and to fermions in Minkowski-space. In the context of FRW 
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cosmology it has been used in a similar way in |3B|. Adding the term Mg on 
both sides of the mode function equation it takes the form 



57^ + " 



fcO 



UJt) 



-V{r)U,iT) 



(5.1) 



with 



V{t) = M\t)-MI 



n 



fcO 



k^ + M^ 



1/2 



(5.2) 



(for the definition of M'^{t) see eq.( |3.17| )). Including the initial conditions 
(|3.22| ) the mode functions satisfy the equivalent integral equation 



f/,(r) = e-*^^«^ + jdT'Ak,rotir - r')F(r')f/fe(r') , 





with 



1 



^fc,retlT - r 



n 



fcO 



■0(r-r') sin (fi,o(r-r')) 



(5.3) 



(5.4) 



We separate Uk{T) into the trivial part corresponding to the case V{t) = 
and a function /ia,(t) which represents the reaction to the potential by making 
the ansatz 

f/fc(r) = e-*^-'^(l + /ifc(r)). (5.5) 

/ife(T) satisfies then the differential equation 



h{r) - 2tQlh{T) = -V{t){1 + hir)) 



(5.6) 



with the initial conditions hk{0) = hk{0) = 0, and the associated integral 
equation 

T 

hkir) = fdr' Ak,rot{r " r')V{r'){l + /ifc(r'))e^^-(--^') . (5.7) 



We expand now hk{T) with respect to orders in V{t) by writing 

h'i\r) + h^P{r), (5.9) 



hk{r) = 4'^(r) + /.f(r) + /.f(r) + --- (5. 
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where hk {t) is of n'th order in V{t) and hj^ (r) is the sum over all orders 
beginning with the n'th one: 

_ °° 

hf\r) = Y.h^^\r). (5.10) 

l=n 

The hf^ are obtained by iterating the integral equation ( |5.7| ). The function 
h^ (t) is identical to the function hk{T) itself which is obtained by solving 
(|5.6| ). The function h^. (r) can again be obtained by iteration via 

T 

hf\r) = /rfr'A,,,et(r - T')V{r')h^P {r')e'''^^'^---'^ , (5.11) 



or by using the equivalent differential equation 

h'fir) - 2tn,ohfiT) = -V{r)hf{r) . (5.12) 

This iteration has the numerical aspect that it avoids computing h\. via the 

small difference h\^ — hk ■ 

The integral equations are used in order to derive the asymptotic be- 
haviour as fifco -^ C)0 and to separate divergent and finite contributions. We 
will give here the relevant leading terms for hi. (r) and h^. (t). We have 

r 

4^)(r) = -^ fdr'{expi2tQ,oir - r')) - l)r(r') . (5.13) 

Integrating by parts we obtain 

/drV(r')-^F(r)+^/rfr'exp(2^fi,o(r-r'))\/' 

(5.14) 



h^k\r) = -^^ J dr'V{r')-^^V{r)+^^ I dr' exv(2zQ,n(r-r'))V'(r' 



or, by another integration by parts, 



h'^\-) = -^ /^^'^(^') - I?^^(^) + ^^'(^) (5.15) 



ifcO ^ ^iifcO ^"fcO 



8fi|o, 



T 

/"dr'exp(2ififco(r - t'))V"{t') . (5.16) 
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We will need often the real part of h\. for which we find 
Here we have introduced the notation Q 

r 

C{f, t) = Jdr' cos(2fifeo(r - r'))/(r') . (5.18) 



which will prove to be useful later. For the leading behaviour of hi. (r) we 

find 

t' 



KU Q Q 



^^ Jdr' JdT"V{T')V{T") + 0{n-,,') . (5.19) 



In terms of this perturbative expansion we can write the mode functions 
appearing in the fluctuation integrals in the equation of motion and in the 
energy-momentum tensor as 

\Uk\^ = 1 + 2Re/if^ + |/ipp , (5.20) 

and 

-tQko(-'^il^K''^^ -2tlmh^P*hfA . (5.21) 
As the potential is real, the leading behaviour of the sums is 

l + 2Re4^) + |/.f)p = l-^V{T) + ^sm{2Q,oT)V'{0) + ^V"{T) 

^ cos{2n,oT)V"{0) + JrrV'ir) + 0(n,t) , (5.22) 



S^feo S^fco 



and 



) 2il,^Q 



- 22lm hl'^ - 2zlm h^'^^h!','^ = —V{t) - ^^ sin(2fi,or)r'(0) - ^^V'iT 



f^^cos(2f]fcor)l^"(0) - ^V'ir) + Oin^,'] 



4fi|o " ^ ^ ^" 4fi3^ 



(5.23) 



^For the numerical computation we use the addition theorem to spht the integral into 
two integrals whose integrands depend on t' only; these can be updated easily. 
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From the Wronskian relation 

UuUl' - U'.Ul = 2«fi,o (5.24) 

we obtain 

2iVLkQ (/IRe h^^ + \hfA - 22lm /if ^ - 2ilm hf*h!^^ = , (5.25) 

which proves to be useful in simplifying the mode integrals occuring in the 
energy- momentum tensor. Using the Wronskian relation we obtain for jf/^P 

\UX = ^l-^l (2Re4'^ + |4'^p) + \K\' . (5.26) 

By means of the equation of motion and Eq. ( |5.26| ) we have 

l^\Uk\" = K\'-m'\U,\' (5.27) 

= - (\/(r) + 2nl) (2Re/ip + |/.«p) + \K\' - Vir) 

6 Renormalization of the equation of motion 

Having expanded the mode functions perturbatively, we are able to sepa- 
rate the divergent parts of the mode sum in an analytic way, leaving the 
finite parts for numerical computation. This allows for a free choice of regu- 
larizations. Furthermore, the analytic expressions for the divergent parts, as 
separated from the mode sum, essentially have the standard form as obtained 
from Feynman graphs, so a comparison with purely analytic approaches is 
straightforward. 

Among the regularizations used in field theory in curved space are: point- 
splitting, dimensional, and adiabatic regularization. Point-splitting regular- 
ization is technically involved as it requires performing the delicate and non- 
covariant limit x' -^ x. Adiabatic regularization actually is a subtraction, it 
has often been used [^, most recently in [|1^]. Adiabatic regularization is 



considered to be well suited for numerical computations, as the entire diver- 
gent part is subtracted, encompassing this way the problem of regularization. 
Here we choose dimensional regularization as it fits in the most appropri- 
ate way into our formalism. We have already used it in our previous work 
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on FRW cosmology. Here we have to correct for some omissions in the fi- 



nite terms, and in particular for the conformal anomaly. In [^ we have 
performed the renormalization by considering the equation of motion for the 
condensate and of the Einstein equations. This is possible in 4-dimensional 
conformally flat space, and leads to a conserved energy-momentum tensor; 



it does not correspond, however, to a properly renormalized local action |30 



and does not work for more general metrics. This shortcoming has been 



criticized in [29 



In our previous publication we have performed the conformal scaling in 
four dimensions, and applied dimensional regularization to the resulting 'flat 
space' equations in analogy to the Minkowski case. We thereby missed terms 
arising from the dimensional continuation of the conformal rescaling, result- 
ing in particular in the absence of some terms involving loga(t). As to the 
equation of motion for the infiaton field, these are the only corrections. The 
renormalization of the energy-momentum tensor is more subtle and will be 
discussed in the next section. 

Using ( p.l7[ ) we split the fluctuation integral of the equation of motion 



into a divergent and a convergent part. Using ( |5.17|) we 

d^'-'k l + 2Reh'i\r) + \h^\r) 



Hr) 



(27r)-i 2(]fco 

d'^~^k 1 / 1 1 

+2Re/^f(r) + |4T)(r)p) . (6.1) 

The first two terms in the integrand have to be regularized, as they lead to 
a logarithmic and a quadratic divergence. We can do another integration by 
parts of the C{V', r) and arrive at 

r rl"-~^k 1 f 1 1 

_ J_ cos(2n„r)V"(0) + ^ - ^<^ + 2ReAf + Ift.l^l 



S^feo S^fco S^fco 



ifc 



Those parts of the integral that involve sin(2r2fcor) and cos(2r2fcoT) develop a 
non- analytic behaviour as r ^ 0. This is due to the fact that the modula- 
tion of the integrand by the trigonometric functions disappears in this limit. 
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Here the non-analyticity is of the form rlogr and r^logr, respectively. In 
the energy- momentum tensor the analogous terms result in 1/r and logr 
singularities. These singularities will be discussed later, in section ^. 

We first rewrite the basic equation of motion, including appropriate counter 
terms, as 



if" + a^ [m^ + Sm + {^-^n + S^ R] ^ 
+ — ^(a/^)V + -^{afiYJ' = . 



(6.3) 



Next we separate from the term XipJ-'/2 the dimensionally regularized diver- 
gent parts. The two relevant expressions behave as 

X r d^-'k 1 AMqV 



(af^y^'^j 



(27r; 



3-e 



2 [k2 + M 



211/2 



327r2 



[Lo + 1) 



and 



-lia,rx.mr)l^^^^^^ 



XipVir) 
327r2 



Lo 



in the limit e 



0. Here we have defined 
Lo 



2 , 47r/i2a2(r) 

7 + ^^^W^-^- 



(6.4) 



(6.5) 



(6.6) 



This cannot be absorbed, e.g., in an MS scheme, by a renormalization 
counter term, as it depends on r, and on the initial mass Mq. Recalling 
that V{t) = M^{t) - Ml and 



M'^ir) 



^'+(e-g)i?(r 



+ 



X{afiY 



n 



12 {n - 1) 



a^R 



we find that the two divergent terms combine into 

-Lo- 



XM'^{t)^ ^ AMqV 



327r2 " 327r2 ' 
The divergences can now be cancelled by the counter terms 

dm — 
5X 
5i 



327r2 
3A2 



L. 

327r2^ ' 

327r2 



(6.7) 

(6.8) 

(6.9) 
(6.10) 
(6.11) 
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where 



2 , 47r/i2 
- + ln— V 



7 • 



e 771^ 

This leaves some finite parts proportional to L — Lq 



In 



mT' 



(6.12) 
and an addi- 



tional finite term due to {n — 4) part of M'^{t) multiplied with the l/{n — 4) 
term of the divergent integral. Finally the renormalized equation of motion 
in n = 4 dimensions reads 



^" + a^ 



m 



Am^ + K + A^ 



1 



R 



V- 



A + AA 
6 



where jFgn is defined as 



Ti 



fin 



a'^R 

'2887r2 



+ ^f 



fin 



with 

^fin 



< 



rf3fc 



1 



-,C{V',t) 



167r2 y (27r)3 2fifco V2^fco 
The finite mass and coupling constant corrections are defined as 



Am^ 
AA 

Ae 



327r2 
3A2 



In 

327r2 
327r2 






In 



rn?a? 
Mf 



<^' + ^'^-^fin = (6.13) 



(6.14) 
(6.15) 

(6.16) 
(6.17) 
(6.18) 



2Re# + |/.?P 



Note that these are time- dependent, due to the occurence of a{T). This is a 
feature we have missed in ||25[ , as we did not continue the conformal rescaling 
to n 7^ 4. 



7 Renormalization of the energy- momentum 
tensor 

In order to derive a renormalized form of the Friedmann equations we have 
to renormalize the energy-momentum tensor. In principle this has been dis- 
cussed long ago and indeed the divergent parts will be found to be in one-to- 
one correspondence to those given in the standard literature [^|. However, 
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we have to discuss this subject in the framework of nonequihbrium quantum 
field theory, and we are interested in particular in the precise form of the 
finite parts which will be the subject of a numerical computation. 

Our previous discussion of the renormalization of the energy-momentum 
tensor had three shortcomings: 

(i) as for the equation of motion we have not extended the conformal 
rescaling to n 7^ 4, these will result in loga(t) contributions, as above ; 

(ii) we have not continued the tensor structures to n 7^ 4; these terms, 
proportional to (n — 4), will be explicitly displayed below, they contribute 
finite terms when multiplied by l/(n — 4); 

(iii) finally, we have discussed renormalization using the Einstein equa- 
tions, and not the Hilbert-Einstein action; this is the reason why the confor- 
mal anomaly did not appear. 

As, nevertheless, the energy-momentum tensor was found to be conserved, 
these shortcomings passed unnoticed. This fact, and the absence of the 
anomaly, are related to the very special nature of the conformally flat metric 
in n = 4. Indeed a unique determination of the counter terms would require 
the consideration of quantum field theory in a g'enera/ background metric, and 
will necessarily be incomplete if only one (and even very particular) metric is 
used. We will first derive the divergent terms and discuss this rather subtle 
matter afterwards. 

In deriving the divergent and finite parts of the energy-momentum tensor 
we have to consider an expansion in powers of (n — 4) = — e. The general 
expression contains various factors with an explicit n-dependence. We define 
a tensor TJ^^f) where in all of these factors n is set equal to 4, in which, however, 
the fluctuation integral is defined in n dimensions. We define a second tensor, 
^h"~^^ which takes into account terms linear in (n — 4) from the explicit 
factors, multiplied by terms proportional to l/(n — 4) from the divergent 
fluctuation integral, and of the counter terms. The remaining terms, less 
singular or with higher powers of (n — 4), do not contribute when n — *> 4. We 
expand the tensor as 

T^. = T^u ('^ + T^u ("-'^ + 0{n - 4) . (7.1) 

Explicitly these tensors are given by 

rr (A) ^ /2 , 1 / 2 , X 2\ 2 , ^ + '^^ 4 
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1 



d^'-'k 1 ( 1 ,,.,,2 , 1 ^, ,,,,,2 1 f, 1\ o,rr,2 






T-T (n— 4) 






2a2 



+5tt I -2^ - SiJ^i?] + 57 

^ \ 2a 72 2 2 / 



2 
i?2 



H^R + 3iJ^ 



+ {n-A)J 



Jn-l 



k 1 1 



(27r)"-i 2fio a 



2-e 



2e 



/J2 + 



i? 



H d 



+ ':^^\Uk? 



H 



R 



(««-') x-s I''' 



2a dr 

The divergent parts of the fluctuation integral are 
1 /• d^'^^k 



S. 



divjfluc 



-.n-2 



(27r)"-i 



2a2 4fifcoa2 ^ ^ 



16fiL«' 



4(6e-l)(f + i/^ 



21] 



fcO 



\/^M 



v^( 



If/A. 



(7.2) 



4fi|o 



r 



4(««-)!4^'''^' 



(7.3) 



Dimensional regularisation of the first three terms in the integral yields 



d''-^k 1 



(27r)"-i a"' 
MHt) 



Vt 



kO 



-Vir) 



le^L 



QAn^a 



2^4 



Lo + 



2 4fi,o 



V'^ir] 



128TT^a 



2^4 



3277^0 



2^4 



m^ + (e - g) ^ + JV^' 
64^^2 



-i^o + 



M4(0) M2M2(r 



128TT^a 



2^4 



327r2a 



2^4 



(7.4) 
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The terms proportional Xirn^ip^ and A^/4(y9^ in ( [7.4| ) are cancelled by the mass 
and coupling constant counter terms. The divergent term which depends on 
m^ determines the cosmological constant counter term, that is 



5k 



m 



QAti' 



:L . 



(7.5) 



The remaining terms in (|7.4| ) combine with the corresponding expressions of 
the remaining parts of i^divfluc: 



--(6e-l)(- 



H' 



-,n-2 



dJ'-^k 



(27r)"-i 



1 



1 



2fi 



fcO 



4fi|o 



V{'i 



m 



6 / 327r2 " 



M 



' 



(7.6) 



and 



-^K 



H 

78— n 



rf"-lfc 



(27r) 



' 41^1.0 



V'(t) 



■32^^^^ 



2^ 



A)aHif^ 



l)-Lo 
a 

n - 



2aHM^ + ( ^ - - ) a'^R' + X^p^p' - XaHip^ 



36 



(a^R' + 2a^H 



+ 0{n - 4) 



(7.7) 



The term A(^ — l/Q)R(p'^ in ( [7.4|) is cancelled by the same term with opposite 
sign in (|7.6|). The H'^ip^-teims in ( [7.6| ) and ( [7. 7) ) are absorbed by the counter 
term proportional to 5^, and the divergence proportional to y^yj' is compen- 
sated by this counter term as well. The remaining (/^-independent but still 
time-dependent divergent terms can absorbed into the counter terms 5a Hu 
and 5Z Gtt- We choose 



5a 
5Z 



(e 



1\2 



327r2 



m 



167r2 



(7.8) 
(7.9) 



This way all divergent integrals appearing in the unrenormalized fluctuation 
integral of the energy are removed by the corresponding counter terms, the 
renormalized expression for the energy will be given below. 
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We shall now comment on the terms leading to the trace anomaly. Ein- 
stein's equation and therefore the energy-momentum tensor contain the terms 



(5a«ff„, + <5^(2)^ ^^-^ 



liv 



(7.10) 



All the three tensors are conserved; furthermore, in 4-dimensional confor- 
mally flat space they are linearly related by Eqs. (|2.11| ) and (|2.12|) , and are 
in fact proportional to each other. Therefore, the condition of a finite energy- 
momentum tensor cannot determine the three coefficients 5a , 513 and ^7, if 
one considers just this restricted class of metrics. So we have to reconsider 
our previous choice of 5a. As we do not consider a more general metric here, 
we have to supply some additional information. 

The divergent part of the energy-momentum tensor for a general metric 
in dimensional regularization has been worked out in [Rl] and is discussed in 



detail in 30 . One finds 



{T,u) 



fMU /div 



327r2 



J_rr 

90 ^'^ 



90 ^^ 



+ u- 



(1)^ 



flU 



(7.11) 



As ^"^'H ^y = H ^i,^ m. n = A and for our metric, our divergent part is con- 
sistent with this expression. It is apparent that within our framework there 
is no way to determine the coefficients of these two tensors, which are the 
source of the conformal anomaly. Therefore, the absence of these terms in 



our previous publication is not, as suggested in [g9|, due to an inconsistency 
of our perturbative expansion. They have to be taken over from a more 
general analysis. Having chosen, accordingly. 



5(3 = -57 



28807r2 



(7.12) 



we continue from ra 7^ 4 to obtain the finite terms arising from an l/{n — 4) 
in the renormalization constants multiplied with the (n — 4) parts of the 
tensors, as displayed explicitly in Eqs. ( p.22[]2.27| ). These finite parts then 
contain the conformal anomaly, and further terms proportional to ^ — 1/6. 
The anomalous parts of the zero-zero component and of the trace are 



5p^'^Hu + 5^Hu 
1 




1 



28807r 



12 



(7.13) 
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and 



rpfi ano 



1 



28807r2 \ a 



Till T3I \ 

+ 2H—-2RH^ + 12H^ 



(7.14) 



After regularization and renormalization the energy-momentum tensor is 
given by 



^ren 



1 /2| 1/2|A 2n2|^ + ^^ 4 

2a4^ 2a2^ ^^ 4!a4 



-V? 



1 



6 



H 



6^-7 + ^0 hr^r-^w' 



2a^ 



+AA + A5 (i)iJii + /\ZGtt + T; 



iq,fin 



(7.15) 



The renormahzed expressions the tensors ^^^Hu and Gu are the ones in four 
dimensions. The finite remnants of the divergent parts have the coefficients 

2 



Aa 

AA 

AZ 



i 



327r2 

4 

m 
647r2 



I) j^m^a^ 



In 



M2 ' 



m 



167r2 



In 






(7.16) 
(7.17) 

(7.18) 



which again are time-dependent, due to the occurence of a{r). 

The fluctuation part of the energy density expressed through flnite mode 
integrals is given by 



T, 



q,fin 



4in,fi„ + nr)^ + i(6e-l)§^L 



M 



R 



H' 



^f 



fin 



T. 



add 



where 

^kiiijfin V j 



3M^ 



1 



d^k 



USn'^a^ 2a* J {2Txf2VLQ 
23 



\h 



(1)'|2 






(7.19) 



(7.20) 



and 



-I add 



ry-iano I 
4 



H 



iaHM'^ - V 



IGtt 



967r2a3 
1 / l\/ Wgf,W ^ 6g^M^ ^ 



or, explicitly, 



T-iadd ry-iano _i_ 

^tt — -'-tt ~^ 



(«-i 



iTC^ 



6 



i?'i/ 



36 



72 



6a 



a" 



3i/2 



(7.21) 



3l^-6j^^-2 7^^^+^^ 



72 



6 



m 



^^-^/f^ (^' - lai/c,) gg^, 



H' 



(7.22) 



Next we have to consider the renormalization of the trace of the energy- 
momentum tensor. We introduce the available counter terms into the un- 
renormalized expression for T/^ and separate the trace according to Eq. ( [7.1D . 
T/^ ^^' is given by 






-,'2 



+2(m^ + <5m2)^+^ + ^^"^ 



(fif 



Qa^ 



-V 



+45A + SZC;^ (^) + (5a ^^^i/,'; ^^^ + T,^ ^°° 
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rf"-iA; 



(27r)«-i 2fi,o 



(6^ 



l?7, 



■/1 2 



-f]^ 



r 



\Uk\ 



+ {h'-R/6)\U,\'--^\U,\ 



.P.^^|..r 



(7.23) 



We can split the trace of the stress tensor into a divergent and into a conver- 
gent part. The divergent part of T/^ ^^^ reads after dimensional regularization 



rpH 

A* div 



327r2a2 ' 
1 



60 



r'(r) 



A 



a 



(7.24) 



m^ + ^^^-(l-60(i/^-f) 



MHT)Ln . 
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Inserting the derivatives of the potential V{t) we see that again the counter 
terms absorb all divergent terms in the fluctuation integral of Tj^. The renor- 
malized trace of the stress tensor takes the final form 



rri^ ren 



-[l-6(e + A0] 



V 



12 



IP 

a 



2 <^' - 2^W 



H 



Lpip +2(m +Am)^-- + 
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+4AA + Aam, + AZG^, + T,^ '^'^'' . 



(7.25) 



The finite fluctuation parts of the trace of the energy-momentum tensor are 
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+ 



2a4 
A_ 
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or, exphcitly, 
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^47r2 l^ 6 



^-r 



a 2 8 . 
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(7.28) 



Comparing to Eq. (3.17) of |^ some differences in the coefficients can be 
absorbed into a different choice of the finite parts of the renormalization 
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constants. The terms proportional to In a(r) given there are included here 
into T^^^^"^ via the coefficients AA, AZ, and Aa; they are identical. We have 
checked, using MAPLE, that the energy momentum tensor is covariantly 
conserved. This is a valuable cross check, as it relates the various terms in 
Ttt and T^^ in a rather complex way. 

8 Removing initial singularities 



The set of initial conditions used in the previous section leads, after renormal- 
ization, to singularities at r = in the remaining fluctuation integrals occur- 
ing in the equation of motion for the infiaton and in the energy-momentum 
tensor. This means that these initial singularities affect the Friedmann equa- 



tions as well. As already shown in [^ these singularities can be removed 
by a suitable Bogoliubov transformation of the naive initial state, i.e., the 
vacuum state of a Fock space based on the mode functions with the initial 
conditions Eq. (|3.22|) . 

Choosing such an initial state one finds, in the fluctuation integral J-gn 
the terms (|6.14|) : 
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sm(2VLkQT) 






cos(2r2fco'7") 



(8.1) 



While these are only nonanalytic as r — i> 0, their ffist and second derivatives 



are singular. They occur in the energy-momentum tensor via ( [7.19|) and 
(|7.26| ). The singular behaviour is given by 
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(8.2) 



(8.3) 



A Bogoliubov-transformed initial state (which is the Heisenberg state of the 
system) can be defined by requiring 



a(k) — pfca^(k) \i) = 
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If the fluctuation integral, the energy and the pressure are computed by 
taking the trace with respect to this state the functions Ukir) get replaced 
by 

Ffc(r) = cosh(7fc)f/fc(r) + e'''^ sinh(7fe)f/*(r) , (8.5) 

where 7^ and 6^ are defined by the relation 

Pk = e'^ tanh(7fe) . 



The fluctuation integral then becomes 
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Expanding as before we find 
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Requiring that the terms proportional to V^'(O) and V'(0) vanish leads to the 
conditions 
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For large k this behaves as 
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The factor cosh(27fe) is equal to 1 for 7^ = 0; the difference w.r.t. 1 is given 
by 

,2/ 



cosh(27.) - 1 = 2sinh^(7.) '^°° ^^ 



(8.12) 
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The new terms proportional to sinh(27fc) behave as behave as 

k^oo \V'{0)\ 
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The dimensionally regularized fluctuation integral (|8.8|) takes, after cancel- 
lation of the singular integrals induced by Eqs. ( p.9| ) and ( |8.1CI| ), the form 
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Using ( ^.12|) and (|8.13| ) we see that the Bogoliubov transform does not affect 
the ultraviolet divergences, and that, therefore, the renormalization proce- 
dure remains unchanged. 

The structure of the equation of motion and the energy-momentum tensor 
remains the same after the Bogoliubov transformation, if the fluctuation 
integral JFgn is replaced by 
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(8.16) 
(8.17) 

(8.18) 



+ sinh(27)ir3(A;,r) + K4{k,T) cosh2(27fc) 
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where K-i{k,T) denotes 

K3(A;,r) = Re{e-2*^'=«^-*'^[/i'^-2^f],o(l + M/i'fe]} , (8.19) 
1 



K,{Kt) 2 



|.(T)'|2_V^^W 



(8.20) 



With these replacements in the equation of motion and in the energy-momen- 
tum tensor the inflaton eqution of motion and the Friedmann equations are 
ultraviolet finite and tree of initial singularities. We should mention that this 
initial state is not determined uniquely. Any state or density matrix based 
on this new 'vacuum state' is admissible as long as the spectrum of coherent 
or incoherent exictations decreases stronger than the one parametrized by 
the Bogoliubov angles •jk and 6k- 

We should like to mention that fixing the initial conditions results in a 
selfconsistency problem, as the finite, modified fluctuation integrals in the 
equation of motion and in the energy-momentum tensor do not vanish at 
r = 0. They depend on the other initial parameters H{0),R{0),ip{0),ip'{0), 
and v'"(0), and vice versa. For the typical parameter sets, this selfconsistency 
problem can be solved by an iteration which converges quickly. 

9 Large N model 

9.1 General formalism 

We now consider the 0{N) a model defined by the Lagrangian 



whith N real scalar fields 0*, i = 1,..,A^. The nonequilibrium state of the 
system is characterized by a classical expectation value which we take in the 
direction of 0Ar. We split the field into its expectation value, or mean field, 
(p and the quantum fluctuations ip via 



0*(x, r) = 5VViV0(r) + ^Xx, r) . (9.2) 

In the large- A^ limit one neglects, in the Lagrangian, all terms which are not of 
order N . In particular terms containing the fluctuation ip^ of the component 
0^ are at most of order yN and are dropped, therefore. This is in contrast 
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to the Hartree approximation where the fluctuations of (p^ are included. 
The fluctuations of the other components are identical, their summation 
produces factors A^— 1 = N{l+0{1/N)). The quantum fluctuations are again 
decomposed into mode functions via ( p.l9| ). The nonequilibrium equations of 
motion for the field 0(r) and of the mode functions Ukir) have been derived 
by various authors [^ ^ . These equations differ from those of the one-loop 
approximation by the fact that the fluctuation integral not only modifies the 
mass of the mean field, but also the one of the quantum fluctuations. 
The equation of motion for the mean field can be written as 



0"(r) + M\t)(P{t) 
the one for the mode functions as 







U'aT)+ e + M\T) U,{r)=0 



(9.3) 



(9.4) 



with the same effective mass 

= a{Tf [m^ + 6m^ + (^ + 5^ - Cn) R{r) 
+ (A + 5A)a(r)^(0(r)2 + ^(r) 



M^(t) 



(9.5) 



We have included here the renormalization counter terms. J-'{t) is again the 
fluctuation integral 

nr)=f^^:^m'. (9.6) 



J (27r)"-i2no 

We restrict ourselves to a temperature T = system, here. Generalization to 
thermal systems is straightforward [^]. As in the one- loop case we rewrite 
the mode equation in the form 






+ n' 



fcO 



Ukir) 



-V{T)Uk{r) , 



introducing, thereby, the time dependent potential 

V{t) = M\t) - M\0) , 



(9.7) 



(9.8) 



The "initial mass" M(0) = mo will be determined below, as a solution of a 
gap equation. We define the time-dependent frequency fikir) via 



^lir) 



e + MHr) 



(9.9) 
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In the large- A^ limit, and in n dimensions the energy density is given by 
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(9.10) 



and the trace of the energy-momentum tensor takes the form 
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(9.11) 



9.2 Renormalization 

The way in which the renormalization counter terms are determined in the 
large- A^ case has been described in detail in [^. We will closley follow this 
approach. As in the one-loop case use the expansion of the mode functions 
in order to single out the divergent parts: 



rrin 



T=-M\T)x,-i^,+T, 



167r2 



fin 



(9.12) 
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with 



Xo 



167r2 



2 , 47r/i2 

7 + ^"m^-^ 



The finite part of the fiuctuation integral reads 



^f 



fin 
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(27r)"-i2fifco V2f] 



^^-C(V',r) + 2Re/if + |4^)|2 



feO 



(9.13) 



(9.14) 



The renormahzation conditions for mass, couphng to the Ricci scalar and 
coupling constant are obtained from the requirement that the frequencies 
which appear in the mode equations are finite, i.e. that M'^{t) is a finite 
quantity: 



m^ + 6m^ + {^ + S^- l)R 



Q' 



+ (A + SX) U^ + T 



m^ + (e - -)R 



+ AU^- 



IGvr^ 
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fin 



(9.15) 



From this condition, and using Eq. (|9.12|) , we find the following counter 
terms 

}?x 

(9.16) 






1-Xx 
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1-Xx 
A 



{x + pX) 



Xx 



e-g)x + (T 



(9.17) 
(9.18) 



The occurence of the free parameters p, a shows that the counter terms are 
not determined uniquely; the choice of these parameters corresponds to the 
freedom of an independent choice of renormahzation conventions. 
We obtain for M^ the manifestly finite expression: 
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Here 
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A 1 m^aJTy ' 
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(9.19) 
(9.20) 



and 
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^.ni^-; 



m^ 



We are now able to determine mg. Taking Eq. 
implicit equation for tjiq = M(0), the gap equation 

1 



(9.21) 
at r = we obtain an 
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With Eqs. (|91^) and (|9]2|) the time dependent potential V{t) = AP{T)~ml 
can be expressed in terms of finite quantities. 

Renormalization of the energy-momentum tensor proceeds analogously. 
We insert the perturbative expansion of the mode functions and determine 
the remaining counter terms so as to render the tensor finite. We require 
them again to be independent of time and of the initial conditions. The 
energy density then is given by 
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explicitly displaying the divergent parts. We have defined 
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After inserting the expressions for the coefficients 6X, 5m^ and 5^ we can now 
fix the coefficients of the higher derivative counter terms, the cosmological 
constant and the wave function renormalization so as to obtain a finite ex- 
pression for the energy density. This is a very tedious exercise; using MAPLE 
we find 
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The coefficients of the finite parts of the counter terms are: 
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AZ = ^,;^(m^ + Ap)(6e-l + 6Aa)ln J/ , (9.30) 



AA = -^(,„^ + Ap)^„!!i!^. (9,31) 

o47r"^ ^ '' rriQ 

The renormalized trace of the energy-inomentuin tensor does not require any 
further counter terms, we find 
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9.3 Removing initial singularities 

We again have to remove the initial singularities. The construction of the 
initial state, and the modified expressions are analogous to those of section 
|[ Of course for the potential V{t) one has to use Eqs. ( |9.8|) , (|9.19|) , and 
( 9.22| ). As in the one- loop case the initial parameters have to be determined 



in a self consistent way, here in addition one has to solve the gap equation. 
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10 Numerical results 

We have numerically implemented the formalism developed in the previous 
sections. We present here some results of the numerical computations. We 
restrict the presentation to results for the 0{N) a model. 

We consider the coupled evolution of the scale parameter via the renor- 
malized Friedmann equations ( p. 15 ), ( 2.14| ) with the energy- momentum ten- 



sor given in Eqs. (|9.28|) and ( |9.32|) , of the inflation fleld via Eq. (|9.3| ), and 



of the quantum modes using Eq. (p.l[). The quantum back reaction on the 
inflaton field and the scale parameter is included from the initial time on, 
this start having made possible by the modification of the initial state. 

We show results for different sets of parameters, choosing values that are 
considered realistic for chaotic infiation, similar to those used in [0. We 
do not consider the evolution during infiation, but only during the reheating 
period. This means that the initial amplitude of the infiaton field 0(0) is 
less or at most of the order of the Planck mass. During infiation the low 
momentum modes increase exponentially, requiring a different numerical or 
semianalytic approach [|19], |10], ^. 

We set A = 10~^^ throughout. This is a value considered to be realistic. 
The inflaton mass is currently estimated to be of the order m ~ 10^^ ■ Mp. 
With such a value fluctuations do not develop, even if v?(0) ^ Mp, and the 
evolution is dominated by the classical amplitude. We therefore consider 
smaller values of the inflaton mass, we also vary the initial fleld amplitude. 
In one group of parameter sets we flx the ratio ip{0)/m and vary the ratio 
m/Mp, as it was done in |12[. In another set we start with ip{0) = 2Mp and 



vary the ratio m/Mp. The parameters of the various sets are displayed in 
Table 1. 

In the flrst group with the parameters sets 1, 2 and 3 we flx ip{0)/'m = 
2 ■ 10'^ and choose m/Mp = 10~^^, 10~®, and 10~^, respectively. This entails 
that the initial amplitudes vary along with the masses. The absolute time 
scale is determined by setting m = 1 in the numerical computations. So 
for all Figures which refer to this data set the abscissa is mr. This also 
determines the units of the various physical quantities. 

For flrst parameter set the initial amplitude is extremely small, and so is 
the back reaction on the scale parameter. We have essentially the situation 
of Minkowski space-time. The scale parameter stays almost constant. The 
evolution of the fleld amplitude (f is shown in Fig. 1. The fluctuations develop 
until jFfjn ~ (f"^ and then have a stationary amplitude, as well-known from 
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studies of the large A^ model [^, ^]. We display, in Fig. 2, the fluctuation 
energy and the total energy, which show a small asymptotic decrease due to 
the expansion. The asymptotic value of the ratio p/£ is found to be 1/3, 
corresponding to a radiation dominated ensemble. 

For the second set the evolution of the scale parameter is plotted in Fig. 
3; it is sizeable as to be expected for an initial amplitude V'(O) = 0.2Mp. 
The classical field amplitude, as displayed in Fig. 4 , starts decreasing after 
mr ~ 5; the asymptotic decrease is stronger than expected for the large- 
time behaviour in Minkowski-space time (r~-^^), and so is due, in part, to 
the expansion. The fluctuations, diplayed in Fig. 5 develop again until the 
asymptotic amplitude is reached, at the same time they are red-shifted. The 
ratio p/ £., shown in Fig. 6 becomes ~ 1/3 (radiation dominated) after the 
fluctuations set in, and decreases to zero asymptotically. This has similarly 
been found in [^ for a similar parameter set. The fluctuation energy be- 
comes of the same order, as the total energy, see Fig. 7 . We finally plot, 
in Fig. 8 , a typical energy distribution of the produced quanta, with the 
familar resonance band at low energy. 

For parameter set 3 the field amplitude starts at 2Mp. The scale param- 
eter develops strongly (see below). The fluctuations hardly develop, they are 
redshifted immediately. The evolution is essentially driven by the classical 
field, with {pjS) ~ 0. We show the classical amplitude in Fig. 9., its decrease 
is mainly due to the expansion. 

These results are qualitatively analogous to those of [|I^], where A was 
chosen 10~^^. The main difference is in the time scale for the built-up of 
fluctuations. In [^ the mass m is of the order vA0(O), while here the latter 



quantity is ten times bigger, and dominates in the estimates for the time 
scales. 

We now consider the second group of parameter sets, for which the initial 
amplitude is flxed, (/^(O) = 2Mp. For the sets 3, 4, and 5 the inflaton mass 
then is chosen as m/Mp = 10"'^, 10"*^, and 10~^^, respectively. Numerically 
we have put Mp = 10^ for all parameter sets. So for all Figures which refer to 
this data set, the abscissa is 10"^ Mpr, and the units for all physical quantities 
are fixed correspondingly. 

The data set 3 has been discussed above, it displays a strong expansion 
and fluctuations play essentially no role. The situation changes with de- 
creasing inflaton mass, which for data set 4 is smaller by a factor of 10. The 
universe then still expands strongly, but at the same time we flnd sizeable 
fluctuations. Still the fluctuation energy only amounts to one part in 10^ of 
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the total field energy. For the much smaller masses value m/Mp = 10^^^ the 
expansion rate only goes down by a factor of about 5, the fluctuation energy 
is of the same order as the total energy. 

We show, in Fig. 10, the evolution of the scale parameter for the data 
sets 3, 4, and 5, all of which show a strong increase. The computation of 
the time evolution for data set 3 was stopped at r ~ 9 as the oscillation 
period in conformal time decreases, requiring very small time steps. The 
evolution becomes uninteresting in the present context, as the fluctuations are 
negligible. In Fig. 11 we display the fluctuation energy for these parameter 
sets. We also plot the total energy for set 5, the total energy for the other 
sets is almost identical, on the logarithmic scale. It is clearly seen that for 
set 3 the fluctuations hardly develop and are redshifted immediately. For set 
4 with a somewhat smaller inflation mass the fluctuations evolve but remain 
on the level of one part in 1000. Finally, for very small masses the fluctuation 
energy increases to the same order of magnitude, roughly 90% of the total 
energy. 

The computations were performed on Pentiumll PC's , the CPU time for 
one of the parameter sets as presented above, was 4 — 8 hours. The typical 
time step, on the scale displayed in the Figures, was 10^'^, with an adaptive 
step-size control. The covariant energy conservation was fulfilled to better 
than 1 part in 10^. 

11 Outlook 

We have presented here the renormalized equations of motion for a self- 
coupled scalar field in a conformally flat FRW universe including the quan- 
tum back reaction in one-loop and large- A^ approximations. We have ap- 
plied this formalism to the post-inflationary preheating period. Our results 
are consistent with those of other authors if we choose the same parameter 
sets. We found that for sufficiently light inflaton masses and initial values in 
the range ip ~ Mp one can have a substantial growth of the energy of the 
fluctuations at the same time as a substantial cosmological expansion. Such 
low values of the mass are not in the range of commonly accepted inflaton 
masses; however, this case could be relevant if the quantum fluctuations are 
those of other flelds. 

We have not considered the inflationary stage itself. In this case, for 
1^ = 0, the low momentum modes grow exponentially due to a term —R/6 in 
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the effective mass, even in the absence of spinodal decomposition. This expo- 
nential growth and the extreme red-shifting pose special numerical problems, 
so that semi-analytical techniques become necessary. The collective evolu- 
tion of the low momentum modes ('zero mode assembly') indeed leads to an 



essential simplification |]10|, ^ 



We have considered here the quantum fluctuations of the inflaton field 
itself. It would be interesting to couple other fields, as e.g. fermion or 
gauge fields to the inflaton. This would allow for a more general choice 
of parameters and could lead to interesting phenomena as e.g. fermionic 
preheating ^2, 



It is known since some time |3^, that gravitational gauge invariance re- 



quires the inclusion of the scalar metric perturbations on the same level as 
the inclusion of the scalar field fluctuations. The resonant growth of such 
metric perturbations in conjunction with those of the scalar field has been 
considered recently ||4^. It would be interesting to extend the formalism 



presented here so as to include metric perturbations. 
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Table Captions 

Table 1: Parameter sets. We display the parameters of the parameter sets 
1 — 5. We only show the various ratios. The numerical values of the dimen- 
sionful quantities m, Mp, and v^(0) in the computer code are fixed by the 
numerical value v^(0) = 2 ■ 10^. 

Figure Captions 

Fig. 1: Classical amplitude v^(t) for parameter set 1. 

Fig. 2: Fluctuation energy (sohd hne) and total energy (dashed hne) for 

parameter set 1. 

Fig. 3: Evolution of the scale parameter a{T) for parameter set 2. 

Fig. 4: Classical amplitude v^(t) for parameter set 2. 

Fig. 5: Fluctuation integral jFgn for parameter set 2. 

Fig. 6: Ratio if pressure to energy density, p/S, for parameter set 2. 

Fig. 7: Fluctuation energy (solid line) and total energy (dashed line) for 

parameter set 2. 

Fig. 8: The energy distribution of the fluctuations at r = 20 for parameter 

2. 

Fig. 9: Classical amplitude v^(t) for parameter set 3. 

Fig. 10: Evolution of the scale factor a{T) for parameter set 3 (dotted line), 

set 4 (solid line), and set 5 (dashed line). 

Fig. 11: Evolution of the fluctuation energy for parameter set 3 (dotted 

line), set 4 (dash-dotted line), and set 5 (solid line). The total energy (dashed 

line) is the one for parameter set 5. 
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set # 


A 


m/Mp 


V9(0)/Mp 


1 
2 
3 

4 
5 


10-12 
10-12 
10-12 
10-12 
10-12 


10-12 
10-8 
10-^ 
10-8 
10-12 


2-10-^ 
2-10-1 

2 

2 

2 


Table 1 
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